Scaling in rupture of polymer chains 



OO 

o 
o 



S. Fugmann and I.M. Sokolov 
Institut fiir Physik, Humboldt- Universitdt Berlin, Newtonstrasse 15, 12489 Berlin, Germany 

We consider the rupture dynamics of a homopolymer chain pulled at one end at a constant loading 
rate r. Compared to single bond breaking, the existence of the chain introduces two new aspects into 
rupture dynamics: the non-Markovian aspect in the barrier crossing and the slow-down of the force 
propagation to the breakable bond. The relative impact of both these processes is investigated, and 
the second one was found to be the most important at moderate loading rates. The most probable 
rupture force is found to decrease with the number of bonds as fmax oc const — (\n{N/r)) ' and 
finally to approach a saturation value independent on A''. All our analytical findings are confirmed 
by extensive numerical simulations. 

PACS numbers: 82.37.-j, 05.40.-a 



00 



o 



X3 
O 

o 



> 

(N 

l> 

O 
OO 

O 



X 



Modern developments of micromanipulation tech- 
niques made possible the experiments on the mechani- 
cal response of single molecules to a well-defined load. 
The typical setup here is the load ramping in which the 
force grows linearly in time until the molecule breaks 
or changes its structure. Examples are single molecule 
pulling experiments [l|, l2, Isl, Ij, la, IgJ and the ones onpro- 
tein unfolding [7|] . The dynamic force spectroscopy [71, l8| 
delivering the complete spectrum of the rupture force ver- 
sus the loading rate offers an extremely powerful tool to 
determine the bonds' strengths and gives deep insights 
into the internal dynamics of the structure under study 

13, [3. 
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From the theoretical point of view the rupture of a 
bond under a constant loading rate can be described 
as a thermally activated escape process from a time- 
dependent potential well [13, El, E, Q [H, [3 ■ Dudko 
et al. [ifll approximated the energy landscape close to 
the barrier up to the third order and predicted that the 
rupture force scales like const -fflnr)^''^, with r being 
the loading rate. It is shown |10l|. that particularly in 



the strong pulling limit [15|, [la] this theory reveals bet 



ter agreement with numerical simulations then a linear 
theory given in [sl, [iJ] • 

In application to polymers this theoretical description 
concentrates on breaking of one bond (presumably the 
weakest one) and fully disregards the chain structure of 
the system. On the other hand, the situation of breaking 
of the more or less homogeneous chain has hardly been 
considered. Thus, Ref. '^X^ discusses breaking of a ring 
of identical bonds, i.e., of a chain with periodic boundary 
conditions. We note that although the experimental real- 
ization of this situation is possible in the setup discussed 
in Ref. [17| , the typical case corresponds to a linear chain 
pulled at one of its ends with another end clamped. In 
what follows we consider this typical situation where the 
force is applied to one of the terminals of the chain while 
the other end is grafted to a surface. 

For sufficiently low loading rates and short chains each 
bond experiences the same force and the rupture can oc- 
cur at an arbitrary bond. For a longer chain or for a high 



loading rate the actual force profile along the chain has 
to be taken into account. When pulling at one terminal 
of the chain the time of the order of the Rouse time tr 
|18J is necessary for the stress to propagate through the 
chain; hence the stress in adjacent bonds differs. If bond 
rupture occurs on a time scale shorter than the Rouse 
time, only a part of the chain is under stress and ac- 
counts for the rupture process. Thus the rupture force 
will be crucially affected by the number of monomers in 
the chain. The aim of this work is to investigate the 
impact of the chain length on the rupture force both nu- 
merically and analytically. Our model will correspond 
to a chain of monomers interacting via the Morse poten- 
tial. Otherwise, the model is identical to the Rouse one: 
we disregard hydrodynamical interactions and describe 
the interaction of the monomers with the heat bath via 
independent white noises. 

The monomer- monomer interaction along the chain is 
modelUed by the Morse potential U{q) given by 



U{q) = £ (1 - ^""^'^' ' 



(1) 



with dissociation energy C /{2a) and stiffness Ca. The 
constant loading enters through an additional time de- 
pendent potential of the form L{q,t) — —qRt with load- 
ing rate R. We note that this potential differs from the 
often used harmonic linker potential [lO|, [l^, [ij, [l9[ , but 
preserves the feature of linearly increasing force within 
time we are focusing on. Its advantage is that the ex- 
trema of the overall potential U{q, t) = U{q) + L{q, t) can 
be obtained analytically. To illustrate the viability of the 
effects let us take the experimentally relevant values of 
parameters - C/{2a) -- 10"^^J with a ~ 10^°m"\ the 
friction 7 — lO^^kg/s, and a loading rate R = 10^'^N/s 
and obtain, that the effects discussed below are impor- 
tant for chains consisting of around 300 monomers, i.e. 
already for relatively short ones. 

To proceed, let us first recall the rupture dynam- 
ics of a single Morse bond under constant loading 
[lOl . [ill [la [ig . [20|. At smaller loads the overall po- 
tential has two extrema, a minimum corresponding to 



a metastable state of the pulled bond, and a maximum 
providing the activation barrier. There exists a critical 
load Fi = F{ti) — C/4 for which the extrema merge at 
qi = ln(2)/ck and disappear. In the purely deterministic 
dynamics the Morse bond breaks exactly at ti = Fi/R. 
Since the system is in contact to a heat bath at temper- 
ature T, its overdamped dynamics is described hy q — 
-C/7(l - exp(-ag)) exp(-ag) + R/-it + ,/2k^Tj^^{t) , 
with ^(i) being Gaussian white noise. We introduce 
c = C/j {[c] ~ nm/fis), r = R/^ {[r] — nm/fxs^) and 
/ = F/7 (in the following / is referred to as force). 
The diffusion coefficient is denoted by Z? = ksT/^ 
{[D] = nm? / ^s). Then, the above given experimentally 
relevant loading rate corresponds to r = 10~'^nm/ pis^ in 
our calculations. 

The probability Wi{t) that the bond remains intact 

can be expressed through the following kinetic equation 

iSEillillOl dWAt)/dt = -k(t)Wi{t), with kit) being 

the Kramers rate |2ll . |22| | . Taking f = rt we can rewrite 

the kinetic equation in the form 



dWijf) 
df 



-k{f)W,{f). 
r 



(2) 



The measured PDF for the rupture forces Pi{f) then is 
P^{f) = -dW,[J)/df. 

Under the assumption that / is close to fi when bond 
rupture occurs, it is usual to expand the potential around 
the inffection point qi up to the third order in devia- 
tions from qi [9|, [lOJ, |ll|]. The Kramers rate becomes 



fc(/)=ca/(47r)v'l-/Mexpi 
Solving ([2|) one derives 



Wi{f) =cxp<^ exp 



-c/(3ai?)(l -///,: 



-i 



(3) 



with V = ca^D/iSn) and w = c/{3aD). In the limit of 
small loading rates the most probable rupture force fmax 
follows the scaling relation [10| 



In (v/r) 



(4) 



Let us now turn to a chain of N bonds. One of its ends 
is fixed and another one is exposed to a monotonously in- 
creasing force with loading rate r = const. To get insight 
into the role of the chain we consider first the situation 
when only one bond is breakable and take this bond to be 
either at the fixed or at the pulled end of the chain. The 
rest of the chain is considered as a harmonic Rouse one. 
The influence of the chain on the breaking properties of 
the bond is twofold: First, due to the coupled dynam- 
ics, the overall noise force acting on the monomer stems 
from the whole rest of the chain and is non-Mar kovian. 
Second, since the force does not propagate immediately 
through the chain, a bond at the grafted end of the chain 
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Figure 1: Most probable rupture force fmax as a function of 
the chain length A''. The chain consist of A'^ — 1 harmonic 
springs and one breakable Morse bond, either situated at the 
fixed wall or at the pulled terminal of the chain. The coupling 
parameter of the harmonic springs c is obtained via the first 
order series expansion of the Morse potential given in Eq. [T] 
i.e., c = ac. Superimposed is the value of fmax obtained 
from a single bond pulling experiment. The loading rate is 
r — 10~^nm/fis^. The remaining parameter values are c — 
3.5nm/fj,s, D = 2 x 10~^nm^ /fj,s, and a = 10/nm. Error 
bars indicate the uncertainty due to the binning of numerical 
data. 



experiences at the beginning the force smaller then the 
one that is applied at the pulled terminal. The rupture of 
a single breakable Morse bond at the fixed wall is affected 
both by the non-Markovian fluctuations and by the force 
profile propagation. In contrast, a bond situated at the 
pulled terminal of the chain feels the instantaneous force, 
and the deviations from the single-bond rupture statis- 
tics are solely due to the non-Markovian character of the 
noise. Fig. [T] shows the most probable rupture force fmax 
as a function of the number of bonds in the chain N for 
both situations. The reference value of the single bond 
rupture setup is also shown. The difference between dif- 
ferent situations is evident. Thus, for a breakable bond 
at a wall fmax lies well above the reference value (dotted 
line), while for the breakable bond at a pulled termi- 
nal fmax lies slightly below the value of the single bond 
rupture measurement when the breakable bond is placed 
at the pulled terminal (dashed line). We conclude that 
non-Markovian fluctuations accelerate the rupture pro- 
cess while the retarded force propagation delays it. Since 
the effect of delay is by an order of magnitude stronger, 
we neglect the impact of the non-Markovian character of 
the fluctuations and neglect correlations introduced by 
noise by assuming that the rupture of different bonds is 
independent. 

We pass to a chain of N breakable Morse bonds. Let 
Wi{fi{t)) be the probability of the bond i to be intact. 
Let r be very small, r <ti 1. Then, we can assume the 
forces acting on each spring along the polymer chain are 



virtually the same. This is expected to be true at least 
for a not too large value of N. The probability that all 
N bonds are intact is then given by Wffif) = Wi{f)^ . 
Then, the probability that a bond breaks in an interval 
[/,/ + df] is PNif) = NWi{f)^-^Pi{f) and is given 
by the same expression as Pi with v changed for Nv. 
Together with ([4]) we derive the following limiting scaling 
relation for the most probable rupture force 



J max — Ji 



In [Nv/r) 



(5) 



The experimentally obtained curves of the most probable 
rupture forces are expected to tend to this scaling relation 
in the limit of vanishing loading rates. 

Let us now turn to higher loading rates. The proba- 
bility that all N bonds in the chain are still intact under 
a pulling force at the terminal of the chain being f{t) is 
then given by WM{fi{t)) = Y{tiWi{h{t)) . Passing to 
the continuum limit {i —f x) wc obtain 



Wn {t) = exp ■ 



\ii {Wi{f{x,t)))dx 



(6) 



Since barrier crossing events are very rare, most of the 
time the motion of the monomers takes place close to the 
quadratic potential minima of bond energies. Therefore 
we derive the force profile by considering a semi-infinite 
chain of harmonic springs pulled at a; = with a force 
f{t) = rt, i.e., by solving the following continuum equa- 
tion for the scalar displacement field q{x, t) 



q{x, t) = caAxq{x, t) + rt5{x) 



(7) 



For not too short chains the impact of the clamped end 
can be neglected. 

Then, the force profile f{x,t) connected to q{x,t) via 
f{x,t) = —cadq/dx is given by 



/(C,t) = /W 



1 - erf 
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(8) 

with ^ — xj \fcai. Linearizing this result close to the 
pulled end. 



and using Eqs. ^ and ^ we get: 
W^^(/(i)) = exp|- "^ 



(9) 



r(^,a(i))-r 



3rw2/3y7(i) 



(10) 



Here I = fi^^cair/r, T{b,z) is the upper incomplete 
gamma function T{b,z) = J^ t^-'^ e-* dt, S{N,f) = 
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Figure 2: Most probable rupture force fmax as a function of 
the loading rate r. The length of the chain is A'^ = 100. The 
remaining parameter values are the same as in Fig. [T] 



2Ny/r/caTTfit)/[fJf{t) - 1] and a{t) = w{l - 
f{t)/fi)^^'^. Finally, the probability density distribution 
reads 



PNifit)) = -WNifit)) 



Vl 



6w2/3/(i)3/2 




F(-,a(i)(l + 5(iV,/))^ 



(11) 



N/l 



,-a{t)il + S{N.f))2 



In Fig. [2] we present the numerically obtained most 
probable rupture force fmax as a function of the loading 
rate r for a chain of length N = 100. Superimposed are 
the prediction of Eq. dH) (dashed line), the prediction of 
Eq. ([5]) (dotted hue) , and the most probable rupture force 
derived from Eq. pT|) . In the limit of small r the most 
probable rupture force tends to the prediction of Eq. ([5]): 
Virtually all bonds account for the rupture process of the 
chain. In the opposite limit of very large loading rates 
only a few bonds contribute to rupture (in the extreme 
limit only the one at the pulled end of the chain) and 
the scaling of fmax is given by Eq. (|4|). The crossover 
behavior is very well reproduced by the theory (Eq. (jlip ). 
Small deviations appear for the intermediate values of 
r where the exact force profile along the chain plays a 
role and the linear approximation in Eq. ^ gets slightly 
inaccurate. Higher order corrections might resolve this 
issue. 

In Fig. [3] we present the numerically obtained most 
probable rupture force fmax as a function of the chain 
length N for r = 10~^n'm/ ^s^. Shown are also the pre- 
diction of Eq. ^ (dashed line), and the most probable 
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Figure 3: Most probable rupture force fmax as a function 
of the chain length A'^ for r = 10~^nm/fj,s^. The remaining 
parameter values are the same as in Fig. [l] 



agreement with the outcome of our numerical simula- 
tions, and for r = 10~''nm//xs^ we get Ng ~ 35 (again in 
agreement with simulations, not shown). 

Let us summarize our findings. Compared to single 
bond breaking, the existence of the chain introduces new 
aspects into rupture dynamics, the most important being 
the delayed stress propagation along the chain. We show 
that the most probable rupture force decreases with the 
length of the chain as fmax oc const — (ln(A^/r))^/^ and 
then saturates at the value depending on the loading rate. 
The results of theoretical considerations are confirmed by 
the ones of numerical simulations. 

The authors thankfully acknowledge valuable discus- 
sions with J. Klafter and M. Urbakh. This research has 
been supported by DFG within the SFB 555 research 
collaboration program. 



rupture force derived from Eq. (fTTj) . For not too large 
systems, the predictions of Eq. (O and of Eq. (fTTj) tend 
to each other and both describe well the results from the 
numerical simulations. Moreover, it is evident that the 
scaling relation Eq. ([5]) gives a lower bound for the most 
probable rupture force. Deviations appear for larger sys- 
tem sizes N, where both the numerical results and the 
prediction of Eq. pT|) saturate. We also state that the 
PDF of the actual rupture forces (not shown) is also ad- 
equately described by Eq. pT|) . 

The overall behavior shown in Fig. [3] corresponds to 
a transition from the small- TV scaling regime to a satu- 
ration after some Ns{r) depending on the loading rate. 
The value of Ng can be estimated as follows. A sin- 
gle bond rupture is governed by Pi{f) and rupture oc- 
curs with highest probability at f{x,t) = fmax with 
fmax given in Eq. (g]). The value of fmax can be 
translated into the most probable time to break via 
tmax = fmax/r. The distribution Pi(/) is strongly 
skewed to the left [10|, |ll|, and its variance is given 



by a? = (2^2/27) />-4/3 (In (t;)) 



"^/^ see Ref. Q. We 
then can assume that the rupture of a bond hardly oc- 
curs if the corresponding force is / < fmax — 2cri . The 
portion of the chain Ng in which the broken bond is 
localized is then determined by the condition that the 
force at the bond number Ng is fmax — 2(Ti at the time 
when the first bond is most probably going to break, i.e. 
f{Ng,tmax) = fmax " 2cri. rcsolviug Eq. ^ for the cor- 
responding value of a; = A^s we get 



Ngir) 



(12) 



For r = IQ^^nm/fis^ we have Ng ~ 100, in qualitative 
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